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Abstract 

We prove the existence of explicit linear multistep methods of any order with positive 
coefficients. Our approach is based on formulating a linear programming problem and estab¬ 
lishing infeasibility of the dual problem. This yields a number of other theoretical advances. 


1 Introduction 

In this work we study numerical methods for the solution of the initial value problem (IVP) 

= /(bw) w(^o) = “0 to < t < r, (1) 

under the assumption that the solution is monotone in time: 

||M(f +/i)|| < ||M(f)II \/h>0. (2) 

Here u : R —>• R'”, and || ■ || is any convex functional. The theory pursued herein is also relevant 
when u satisfies a contractivity or positivity condition. We will usually write f{u) instead of 
f{t,u) merely to keep the notation simpler. 

We focus on the class of methods for solving ([T]) that are known as linear multistep methods 
(LMMs). To solve ([T]) by a linear multistep method, we define a sequence of times to < ti < • • • < 
fjv = T where t„ := t„_i + h, and compute the values u„ « u{t„) sequentially by 

k k 

U-n = ^j^n-j A h YYj (3) 

j=l 
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Some prescription must be given for the starting values uq,. . If j6o = 0, the method is 

said to be explicit; otherwise it is implicit. 

We are interested in methods that preserve a discrete version of the monotonicity condition 
©; namely 


||w„|| < max{||M„_i||,...,||u„_fc||}. (4) 

The backward Euler method achieves (IDl under any step size, as long as / is such that (l2)l is 
satisfied. 

In order to achieve the discrete monotonicity property dD with explicit methods, or with 
implicit methods of order higher than one, we assume a stronger condition than (|2)l. We require 
that / be monotone under a forward Euler step, with some restriction on the step size: 

||m +/2/(i/)II < ||m|| VO < ?! </2fe(w). (5) 

Under assumption © it can be shown that any method dUl with non-negative coefficients pre¬ 
serves discrete monotonicity dD under some time step restriction. The step size restriction that 
guarantees monotonicity is then 


0 < Iz ^ C/lpF, 


where the factor 


C ■- 


min a,/ B; 

i<j<k ' ’ 

< 


0 

\ 


if CLj > 0 for 1 < i <k and fij > 0 for 0 < j < k, 


otherwise 


( 6 ) 


is known as the threshold factor or strong stability preserving coefficient (SSP coefficient) of the 
method. 

Runge-Kutta methods with positive SSP coefficient are subject to restrictive order barriers: 
explicit methods have order at most four and implicit methods have order at most six. SSP linear 
multistep methods of high order are therefore of particular interest. The objective of the present 
work is to investigate existence of methods with C > 0 and bounds on the value of C for methods 
with a prescribed order of accuracy and number of steps. 

In this paper we have assumed a fixed step size; the case of SSP explicit linear multistep 
methods with variable step size is covered in IIHKLNI . 


1.1 Previous work 

Contractive linear multistep methods were studied by Sand, who constructed a family of implicit 
methods with arbitrarily high order USanShl Theorem 2.3] whose number of steps is exponential 
in the order of accuracy. Later Lenferink deduced many properties of the optimal methods 
and threshold factors for low order methods HLen89irLen91l . More recently, such methods have 
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been studied by Hundsdorfer & Ruuth ||HR05[ IRH05I , who consider the effect of special starting 
procedures and of requiring only boundedness rather than strict contractivity or monotonicity 
A fast algorithm for computing optimal methods, along with extensive results, was given in 
IIKet09ll . 

1.2 Scope and main results 

Linear multistep methods are closely related to polynomial interpolation formulas. We investi¬ 
gate their properties herein using the framework of linear programming and results on polyno¬ 
mial interpolants. 

The main results proved in the present work are: 

1. Existence of arbitrary order explicit LMMs with C > 0; 

2. A sharper upper bound on C for implicit LMMs; 

3. Behavior of the optimal value of C for k-step methods as k —> oo. 

Along the way, we also give a new proof of the known upper bound on C for explicit LMMs, and 
a new relation between C for certain implicit and explicit classes of LMMs. 

2 Formulation as a linear programming feasibility problem 

By replacing ii„ with the exact solution w(f„) in (|3)l and expanding terms in Taylor series about 
tn, we obtain the following conditions for method (H) to be consistent of order p: 

k 

= l (7a) 

7=1 

^ = 0 1 < ra < p; (7b) 

7=1 ^ ^ 

here and in the rest of the paper, we assign 0^ the value 1. 

A linear multistep method (O has SSL coefficient at least equal to r > 0 iff 

> 0 and a.j — rfij > 0 for 1 < / < k, and jSo > 0. 

Thus the problem of whether there exists a method of order p with k steps and SSL coefficient at 
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( 8 ) 


least r > 0 can be formulated for explicit methods as IIKet09[ LP 2]: 

Find 

>0 Sj>0 

such that 

+>■/!/))"-ft™/"-') =0" 

/-I 

and for implicit methods as IIKet09[ LP 4]: 

Find (9) 

/lo > 0 and j6y >0 Sj >0 for 1 < j <k 

such that 

Y^{Sj + r^j) = l and 

^ ( (Sj + - j6oO'"-i =0 for 1 < m < p. 

j=i 

In both cases, ay can be obtained for a fixed r via ay = Sj + r/ly for all l<i<k. 

3 Optimality conditions for SSP methods 

In this section we develop the basic tools used in this paper. Our analysis relies on Parkas' lemma 
and on the Duality and Complementary slackness theorems in linear programming (LP), which 
we now recall; see, e.g. IISch98l . 

Proposition 3.1 (Parkas' lemma). Let A G Jje a matrix and let b G R"' be a vector. The system 

Ax = b,x >0 is feasible if and only if the system A^y > 0, b^y < 0 is infeasible. 

Proposition 3.2 (Duality theorem in linear programming). Let A G be a matrix and let 

b G R'",c G R” be vectors. Consider the primal-dual pair ofLP problems 

Maximize c^ x subject to Ax = b and x > 0, (10a) 

Minimize b^y subject to A^y > c. (10b) 

If both problems are feasible, then the two optima are equal. 

Proposition 3.3 (Complementary slackness theorem). Let A G be a matrix and let b G R"*, c G 

R” be vectors. Consider the primal-dual pair ofLP problems (llOall and dlObl) . Assume that both optima are 
finite; let x be an optimal solution to (llOall and let y be an optimal solution to dlObll . Then x^ {b — A^y) = 
0, i.e. if a component of x is positive, then the corresponding inequality in A^y > c is satisfied by y with 
equality. 


for 1 < j <k 
for 0 < m < p; 
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Our analysis also relies on the following well-known result of Linear Algebra that can be 
found in many textbooks, for instance, in IISIF08I . 

Proposition 3.4 (Rouche-Capelli theorem). Let A G be a matrix and let b G R"* be a vector. 

The system Ax = b has a solution if and only if the rank of its coefficient matrix A is equal to the rank of 
its augmented matrix, {A\b). If a solution exists, then the set of solutions forms an affine subspace o/R” 
of dimension {n — rank A). 

We can now state three simple lemmas. Since the system (|lla|l - (lllcb below is the Farkas-dual 
of (H), and the system formed by (lllab - dllcll and (llldb is the Farkas-dual of ®, Lemma l3dl and 
Lemma 13.21 are consequences of Farkas's lemma (Proposition 13.1b and the observation that the 
dual variables can be viewed as the coefficients of a polynomial of degree at most p. 

Lemma 3.1. Let r G R and let k, p be positive integers. The following statements are equivalent: 

(i) The LP feasibility problem (|8ll is infeasible. 

(ii) There exist ym G R/or 0 < m < p that satisfy the system 


P 


E ^ 0 

m=0 

for 1 < i < k 

(11a) 

E (-rnr^^ym+rf'ym) >0 

for 1 < i <k 

(11b) 

yo < 0. 


(11c) 

(Hi) There exists a real polynomial q of degree at most p that satisfies the conditions 


‘lii) > 0 

for 1 < i < k 

(12a) 

1 

IV 

o 

for 1 < j < k 

(12b) 

q{0) < 0. 


(12c) 


Lemma 3.2. Let r G R and let k, p be positive integers. The following statements are equivalent: 

(i) The LP feasibility problem ® is infeasible. 

(ii) There exist y^ G R/or 0 < m < p that satisfy the system defined by (lllab - dllcb and 

-yi>o. (lid) 

(in) There exists a real polynomial q of degree at most p that satisfies the system defined by dl2ab - dl2cb 
and 

- q\0) > 0. (12d) 
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Recall that due to the Fundamental theorem of algebra, a real polynomial can be factored 
in a unique way as q{x) = c ■ ~ ^m), where c G R \ {0} is the leading coefficient of 

q and Ai, A 2 ,..., Ap^ G C are the roots of q. Moreover, in view of the Complex conjugate root 
theorem, the non-real roots of q occur in complex conjugate pairs. The proof of Lemma 13.31 is 
also based on the formula logarithmic derivative of the polynomial 

q{x) := A^). 

Lemma 3.3. Let r G R, and let q be a non-zero real polynomial of degree po > 0. Let c G R \ {0} 
be the leading coefficient of q and Ai,A 2 ,.. .,Apg G C the roots of q. For any real y, let s{}i) := 
\{m : 1 < m < po, and Am G R, and Am > }i} \ denote the number of real roots of q that are greater than 
or equal to p. Then the following statements hold. 

(i) The polynomial q satisfies the inequality (|12a|l for some integer j if and only if either of the following 
is true: 


• i is a root of q; 

• i is not a root of q, and 

c. (-i)sO') > 0. (13a) 

(ii) Suppose that q satisfies the inequality (I12all for some integer j. Then q satisfies the inequality (Il2b|) 
for the same value of j if and only if either of the following is true: 


• i is a multiple root of q; 

• i is a simple root of q, and (I13a|) holds for j; 

• j is not a root of q, and 

Po 1 

E ^^ ^ 

(Hi) The polynomial q satisfies (I12cll if and only ifO is not a root of q and 

c. (-i)s(o) < 0. 


(13b) 


(13c) 


(iv) Suppose that q satisfies (ll2cD . Then q satisfies (Il2db if and only if 



(13d) 


For convenience, we introduce the following notation. Let CexpiKp) denote the "optimal" 
SSP coefficient for explicit LM methods with k steps and order of accuracy p, i.e. 

^exp(^/F) — sup C, 

/c-step explicit 
methods of order p 
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and let Ciinp(fc/ p) denote the "optimal" SSP coefficient for implicit LM methods with k steps and 
order of accuracy p, i.e. 

Cimp{k,p)= sup C. 

/c-step implicit 
methods of order p 

Our purpose is to study the behavior of C^xpiK p) and Cimp{k, p). It turns out to be convenient 
to study polynomials that satisfy the following conditions, in place of (I12a|l - (ll2db : 


q[j) >0 for 1 < ; < /c (14a) 

^ ■ qii) — ^ for l < ; < Ic (14b) 

q{0) = 0 (14c) 

q'{0) = 0. (14d) 


Before we proceed to our main lemmas. Lemma 13.61 and Lemma 13.71 we need to prove two 
auxiliary results relating these sets of conditions. 

Lemma 3.4. Let r G K and let k, p he positive integers. Suppose that q is a real polynomial of degree 
0 < po < p that satisfies the system (I12all - (ll2cll . Then there exists a non-zero real polynomial 


V 

q{x) := c • P]^(x - A,„) 

m=l 

of degree p that satisfies the system (I14a|l - (ll4cll . has a leading coefficient c = ±1, and whose roots satisfy 
G "H := |z G C : I Im(z) | < and 0 < Re(z) < lc| for all 1 < m < p. 

Proof Let r G IR, and let k, p be positive integers. Suppose that q is a real polynomial of degree 
Po < P that satisfies the system (I12all - (ll2cl) . Let us consider the polynomial q{x) := q{x) — q{0). 
This is a non-zero real polynomial of degree po that satisfies the system (|14a|l - (ll4cb : furthermore, 
it fulfills the condition q{j) > 0 for all l<i<k. Let q{x) := c ■ 11^=1 factored 

form of q, let M := \{Xm '■ 'i < m < Pq, and A^ G R, and A^ < 0}| be the number of negative 
real roots of q, and let us take 


c := 


Am . — 


q{x) := 


= f — l'jM+p-po . _E_ 

' ' Irl 


if 1 < ra < Po and A„, G TL 
if 1 < ra < Po and A„, ^ TL, or if po < m < 

P 

C- n(x - Am). 


It is clear that q satisfies condition (I14cb , and that c = ±1. 
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Let s and s be defined as in Lemma 13.31 with respect to q and q, respectively. Because of 
statements (E)) and (Iii]l of Lemma l33l we have c • (—> 0 and ~ ^ for all 

l<i<k. Since the non-real roots of q come in complex conjugate pairs, s{j) = s{j) + M + p — po 
mod 2 for all 1 < / < fc. Hence c ■ (—= c/|c| • (—> 0 for all 1 < j < k, and thus 
statement (E)) of Lemma [3.31 implies that q satisfies condition (I14a|) . 

Now we prove that q satisfies condition (Il4b|l . To achieve this, we show that the contribution 
of the roots of q to the left hand side of (Il3bll is less than or equal to that of the roots of q. 

Let us consider first a p E C such that | Im(f/) | > Then 

1 1 _ 2(;-Re(;^)) 1 2 2 

j-p^j-p (/-Re(f/))2 + (Im(/^))2 - |Im(/^)| ^ k-l-j-k 

for all 1 < / < ?c — 1. 

Let us consider now a p E C such that either Re(^) > k or Re(p) <0. Then 

1 I 1 _ 2(;-Re(;^)) ^ ^ ^ ^ 

j-H j-H (; - Re(/i))2 + (Im(^/))^ “ Re(f/) -; j-k 

for all 1 < ; < ?c — 1. 

It follows that 

Pi Po 1 Po 1 

m=l ] 1 ,„=i ] Am 

for all 1 < i < k — 1. If fc is not a root of q then the same holds for j = k; otherwise fc is a root of 
q and c • (—1)®^^) > 0. In both cases, statement © of Lemma 1331 implies that q satisfies condition 
(I14bb . which completes the proof. □ 

Lemma 3.5. Let r G R and let k, p he positive integers. Suppose that q is a real polynomial of degree 
0 < po < P lhat satisfies the system formed by (|12a|l - (ll2cl) and (Il2d|l . Then there exists a non-zero real 
polynomial 

P 

q{x) := c ■ - Am) 

m=l 

of degree p that satisfies the system (I14al) - (ll4dll . has a leading coefficient c = ±1, and whose roots satisfy 
Am G "H := |z G C : I Im(z) | < and 0 < Re(z) < A:| for all 1 < m < p. 

Proof Let r G R, and let k, p be positive integers. Suppose that is a real polynomial of degree 
Po < p that satisfies the system formed by (I12al) - (ll2cl) and (Il2dl) . Let us consider the polynomial 
q{x) := q{x) — q{0). This is a non-zero real polynomial of degree po that satisfies the system 
formed bv (I14al)-(ll4cl) and (Il2db. furthermore it fulfills the condition q(i) >0 for all l<i<k. 

Let q{x) := c ■ Ylm^ri^ ~ Am) be the factored form of q, and let s be as in Lemma 1331 with 
respect to q. Statement © of Lemma |33I implies that c ■ ( —1)®(^) > 0 and statement © of Lemma 
13.31 implies that either c ■ (—1)®^^) > 0 or 0 is root of q of multiplicity at least two. In both 
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cases there are two distinct indices 1 < mo ^ nii < po, such that A,„g = 0, and E R, and 
0 < A„,j < 1. Let M be the number of negative real roots of q, as in the proof of Lemma l3!4l and 
let us take 


c := (-l)^+P^Po . 


Am; !- < 


Am if 1 < m < Po and Am G "H and m ^ mi 
0 if m = mi 

k if 1 < m < Po and Am ^ T-L, or if po < m < p 


q{x) := c • - Am). 


m—l 


We can proceed in the same way as in Lemma 13.41 to show that this q satisfies the claims of 
the lemma. □ 


Lemma 3.6. Let r > 0 and let k, p he positive integers. Then the following statements are equivalent. 

(i) Cexp(k, p) ^ r. 

(ii) There exists a non-zero real polynomial q of degree at most p that satisfies the system (|14a|l - (ll4cb . 

Proof. We can check that the polynomial q{x) := x — ^ of degree one satisfies the system (I12al) - 
(I12cl) with r = 2 for any value of k, and thus Lemma IAT] implies that CgxpiK p) < 2, and hence it is 
finite, for all positive integers k and p. Now we are ready to show the equivalence of statements 
(0 and dull. 

Proof that statement (Hil implies statement 0. Suppose that q := c ■ nm=i(^ “ Am) is a non¬ 
zero real polynomial of degree 0 < po < p that satisfies the system (IMaMMcIl with r = po > 0 
for some positive integers k and p, but Cexp(k, p) = pi > po- Since 0 is a root of q, we can 
assume without the loss of generality that Ai,..., Ap^ 0 and Ap^+i,..., Ap^ = 0 for some integer 
0 < Pi < Po. 

Let us choose 0 < e < 1 so that < PWPO 

, and let us take the polynomial 

L ■ nS.i(^ if c ■ nLiC-.'™) < a 

Let s and s be as in Lemma 13.31 with respect to q and q, respectively. By construction, both the 
leading coefficients and the set of real roots greater than e of q and q coincide, and thus s(/) = s(/) 
for all 1 < i < k. 
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We can check that q satisfies conditions (I12all and (I12cll . Moreover, in view of Lemma [3.31 if j 
is not a root of q for some l<i<k then (I13bll holds for the roots of q with r = po, and thus 


Po 


> 


y 1 


Pi 

E 

m 


1 1 1 

^ + (po - Pi )7 > E 73T- 

=ll ] A„ 


1 

i 


y ^ _l_ _J_ _ ^ > y ^ ^ _ Pi _ 

j-e Hi -e) - j-e 2 


Po + Pi > y ^ ^ 

2 i-e' 


Hence, in view of statement (|ii| of Lemma l33l q satisfies condition (I12bb with r = . Now 

Lemma 13.11 implies that C < for any k-step explicit LM methods with order of accuracy 

p, which leads to the contradiction CexpiKp) < This completes the first part of the 

proof. 

Proof that statement (01 implies statement (El). Let k, p be positive integers. We can check 
that if a non-zero real polynomial q satisfies the system (I14all - (ll4cll with r = p then it satisfies the 
same system with any r > p. Hence it is enough to prove that there is a non-zero real polynomial 
q of degree at most p that satisfies the system (IMaMMcIl with r = C^xpiK p)- 

Let be as in Lemma |3H] By the definition of Cexp(fc/ p) and due to Lemma l3dl there is a real 
sequence and a corresponding sequence of real polynomials of degree at most p, 

such that lim„^ooN = C^xpiK p), and for all n, the polynomial qn satisfies the system (I12all - (ll2cll 
with r = Yn- Due to Lemma l3!4l for any positive integer n, there exists a non-zero real polynomial 
qn{x) := Cn ■ ~ ^n,m) of degree p, such that qn satisfies the system (I14all - (ll4cll with r = r,„ 

moreover c„ = ±1 and G T-L for all 1 < m < p. 

Since for all n, and since W C CP is closed and boimded, the 

Bolzano-Weiertrass theorem implies the existence of a strictly increasing sequence of indices 
such that for all 1 < m < p, the limit A^ := lim£_^oo ^ni,m exists, A^ G TL and non-real A 
occur in complex conjugate pairs, moreover c„i = c for all £. A continuity argument yields that 
q{x) := c ■ ~ Am) is a real polynomial of degree p that satisfies the system (|14a|l - (ll4cb 

with r = Cexp{k, p), which completes the proof. □ 

Lemma 3.7. Let r > 0 and let k, p be positive integers. Then the follozving statements are equivalent. 

(i) Cinip(k, p) < r. 

(ii) There exists a non-zero real polynomial q of degree at most p that satisfies the system (|14a|l - (ll4d)l . 

Proof. The proof is essentially the same as for Lemma 13.61 □ 

The results of this paper are based on analysing the polynomials q appearing in Lemma 13.61 
and in Lemma 13.71 and for this reason, we establish some basic properties of these polynomials 
in Lemma [3.81 and in Lemma [3.91 
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Lemma 3.8. Let k,p be positive integers and let r > C^xpiKp) be a real number. Then there exists a 
non-zero real polynomial q that solves the system (I14all - (ll4cll and has the following properties. 


1. The degree of q is exactly p. 

2. The real parts of all roots of q except for 0 lie in the interval [1, k]. 

3. All real roots of q are integers. 

4. The multiplicity of the root 0 is 1. For 1 < j < k — 1, if j is a root of q then its multiplicity is 2. 

5. For even p, k is a root of q of multiplicity 1. For odd p, ifk is a root of q then its multiplicity is 2. 

6. The polynomial q is non-negative on the interval [0, k]. 

Moreover, ifCgxpiK p) > 0 and r = Cgxp{k, p) then every non-zero real polynomial q of degree at most p 
that solves the system (ll4all - (ll4cD has these properties. 

Proof. Case r = Cgxp{k, p) > 0. Suppose that positive integers k, p are such that Cgxpik, p) > 0. We 
know that Cexp (fc/p) < cxj, see the proof of Lemma l3.6[ Letq{x) := c ■ — ikm) a non-zero 

real polynomial of degree po < p that satisfies the system (I14all - (ll4cll with r = po := Cexp(k, p); 
the existence of such a polynomial is guaranteed by Lemma 13.61 Due to condition (ll4cL = 0 
for some 1 < m^ < po- Now we show that q has all the required properties. 

Proofs for properties |l] and |2j Suppose that either property [T] or |2] does not hold for q. Let 
M := \{m : 1 < m < Po, and m mo, and A„, ^ be the number of roots of q with real 

parts outside the interval [l,k], not counting Amg- Let us take 


c := {-i)M+p-po.c 

{ Am if 1 < m < Po and Re(Am) G [1/^]/ or if m = mo 
k if 1 < m < Po and m f mo and Re(Am) ^ [1/^]/ or if po < m < p 


q{x) := c • J^(x - A„ 


m—l 


By following the argument in the proof of Lemma 13.41 we can show that q satisfies the system 
(I14all - (ll4cll with r = pi for some 0 < pi < po- Now Lemma l3!6l yields Cgxp{k, p) < pi < po/ which 
contradicts our choice of q. 

Proof for property |3j Suppose that property |3] does not hold for q. Let us take 


A 


m 


q{x) 


{ Kn if 1 < WZ < Po Am IR 

[Am\ if 1 < wz < Po and A^ G IR 

Po 

c • rT(x - Am). 

m=l 
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Let s and § be as in Lemma 133] with respect to q and q, respectively. We can check that §{]) = s{j) 
for all 0 < j < k. Moreover, if j is a root of q for some integer 1 < j < k then j is a root of q too, 
with at least the same multiplicity. If, however, j is not a root of q for some integer 0<i<k then, 
assuming that at least one root has moved, we have ~ In view 

of Lemma 13.31 the polynomial q satisfies the system (I14all - (ll4cll with r = pi for some 0 < pi < po, 
which again leads to a contradiction. 

Proof for property m We already know that properties UHl] hold for q. Suppose that = 
A ,„2 = A „,3 = i for some integer l<i<k and indices 1 < nii < m 2 < < p. Then we can take 


Am ■ — 


Am 

Am 1 


if m 7 ^ ra 3 
if m = m 3 



m=l 


A-m)/ 


and, as before, we can check that the polynomial q satisfies the system (I14at - (ll4cll with r = pi 
for some 0 < pi < po, which leads to a contradiction. Hence if j is a root of q for some integer 
l<i<k then its multiplicity is at most 2. 

Suppose now that jo is a simple root of q for some integer 1 < jo < K and suppose further 
that i is not a double root of q for all integers jo < j < k. Let s be as in Lemma 13.31 and let 
jo < ji ^k be the smallest integer that is not a double root of q. Then c • (—l)®*'h) = i)s(yi)^ 

and in view of statements (H) and (Iii]l of Lemma 13.31 q does not satisfy the inequalities (114a ll and 
(Il4b|l for both j = jo and j = ji. Therefore if jo is a simple root of q for some integer 1 < jo < k 
then j is a double root of q for all jo < j < k, and thus Amj = k for some index 1 < mi < p. Now 
we can take 


Am. . — 



if m 7 ^ mi 
if ra = rai 


q{x) := c • J^(x - Am), 

m=l 


and check that q satisfies the system (I14a|l - (ll4cb with r = pi for some 0 < pi < po, which again 
leads to a contradiction. 

Proofs for properties [5]and|6]Property|5]is a consequence of properties [Ull] and lU and the fact 
that non-real roots come in complex conjugate pairs. Property | 6 ] is a consequence of properties |3] 
andlH 

Case r > Cexp{k, p) > 0. We can see that if C^xpiK p) > 0 and q satisfies the system (I14all - (ll4cll 
with r = CexpiK p) then q satisfies the same system with any r > Cexp(k, p), and thus the claim of 
the lemma holds true in this case. 

Case Cexp (k, p) =0 and p < 2k -|-1. In this case, we can construct a non-zero real polynomial q 
of degree p that satisfies the system (|14a|l - (ll4cb with r = 0 and has propertiesHHH by applying the 
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transformations presented in the proof of the first case to an arbitrary non-zero real polynomial 
q of degree po < P that satisfies the same system with r = 0 . 

Case CexpiK p) = 0 and p > 2k + 1. For odd p we define q{x) := x{x — k + (x — k — 

■ Y^j^i{x — iY, and for even p we define q{x) := x{x — k + i)^(x — k — i)^^{k — x) ■ 
~ iY- We can check that these polynomials satisfy all the required conditions. □ 

Lemma 3.9. Let k and p > 2 be positive integers and let r > Cimp{k, p) be a real number. Then there 
exists a non-zero real polynomial q that solves the system (|14a|l - (ll4d)l . and has the following properties. 

1. The degree of q is exactly p. 

2. The real parts of all roots of q except for 0 lie in the interval [1, k]. 

3. All real roots of q are integers. 

4. The multiplicity of the root 0 is 2. For 1 <}< k — 1, if j is a root of q then its multiplicity is 2. 

5. For odd p, k is a root of q of multiplicity 1. For even p, ifk is a root of q then its multiplicity is 2. 

6. The polynomial q is non-negative on the interval [0,k]. 

Moreover, if Cimp{k, p) > 0 and r = CimpiK p) then every non-zero real polynomial q of degree 
at most p that solves the system (I14a|l - (|14d|l has these properties. 

Proof. The proof is essentially the same as for Lemma 13.81 □ 

The following two lemmas provide a way to analyze both the coefficient structure and the 
uniqueness of "optimal" SSP LM methods with a given number of steps and order of consistency. 

Lemma 3.10. Let positive integers k,p be such that CexpiK p) > 0. Assume that and 

are the coefficients of an explicit LM method of order p with k steps and with SSP coefficient CexpiK p), 
and let q be any non-zero polynomial of degree at most p that satisfies the system (|14a|l - (ll4c|) with r = 
CexpiK p). Let us define the sets T := {j : 1 < j < k, (ind qij) = 0} and J := {j: 1 < j < 

K and — q'ii) + Cexpik, p)qij) = 0}, and let us introduce the vectors aj := if^Y^^^for 0 < j <K and 

bj := i—mf^^^ -\- Cexpik, p)fKm=of^^ 1 < / < k. Then the following statements hold. 

(i) If 5j 7 ^ Q for some I <] <K then qij) = 0 . 

(ii) If f>j f D for some I < j <k, then —q'ij) -F CexpiK p)‘lij) = 0. 

(in) If the vectors ajfor ; G X and bjfor j F J are linearly independent then there exists no other explicit 
LM method of order p with k steps and SSP coefficient C^xpiK p)- 

(iv) The polynomial q can be chosen so that the vectors ajfor / G X and bjfor j E J span a p dimensional 
subspace o/RP+^. 
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Proof. The proofs of statements (H) and (HD) of the Lemma. Suppose that positive integers k and 
p, coefficients and and the polynomial q satisfy the conditions of the lemma. Let 

us consider the following primal-dual pair of LP problems. 

Maximize 0 with respect to subject to (jS]) with r = Ceyq,{k, p). (15) 

Minimize yo with respect to {ym)m=o subject to (IllaD and (lllbb with r = C^xpik, p). (16) 

On the one hand is an optimal solution to (ITSll , and on the other hand y^ := 0 for 

0 < m < p is a feasible solution to dlhll . Therefore both (ITSl) and (O are feasible, and thus the 
Duality theorem in LP (Proposition 13.21) implies that the two optima are equal, hence a feasible 
solution (ym)m=o is optimal if and only if yo = 0. 

Let ym := q^"^f0)/ni\ for 0 < ni < p be the coefficients of q. Then y is an optimal solution 
to (O, because = ^(/) ^ 0 fo’’ ^ < j <K and + Cexp(^,p);™)ym = 

-fii) + Cexp{k, p)q{i) > 0 for all 1 <] <k, and yo = q{0) = 0. 

Suppose that 3j 0 for some 1 < i <k. Then the Complementary slackness theorem 
(Proposition 13.31) implies that q{j) = ~ Similarly, if 0 for some \<]<k then 

the Complementary slackness theorem implies —f{j) + Cexp{k,p)q{j) = 0. Hence the proofs of 
statements (HJ and (ju)) are complete. 

The proof of statement duU of the Lemma. Suppose that positive integers k and p, coefficients 
and (^y)y^^, and the polynomial q satisfy the conditions of the lemma. Suppose further 
that the vectors Uj for j E 1 and bj for j ^ J are linearly independent. Consider the system of 
linear equations obtained from the equality constraints of the LP problem dl^ by fixing 5j for 
y 0 X and fj for j ^ J io zero. In view of statements (HI and (ju) of the lemma, this system is 
satisfied by the possibly non-zero coefficients and of any optimal method. Since 

the matrix of this system is formed by the column vectors aj for y G X and bJ for j G JT, 
its columns are linearly independent, and thus the Rouche-Capelli theorem (Proposition 13.41) 
implies that {^j)jej is a unique solution to this system. From this, the third statement of 

the lemma follows. 

The proof of statement (Hil of the Lemma. Suppose that positive integers k and p, coefficients 
(<5;)y^i and and the polynomial q satisfy the conditions of the lemma. Suppose further 

that the polynomial q is chosen so that it satisfies the maximum number of inequalities in (|14a| - 
(Il4bb with equality. Let y,„ for 1 < m < p be as before. 

We show that y satisfies at least one inequality in (|lla|) - (lllb| with strict inequality. Suppose 
that is not the case, then for all l<i<k integers j is a multiple root of q, and thus q satisfies the 
system (I14ab - (ll4cb with r = 0. Lemma 1331 leads to the contradiction C^xpiK p) = 0. 

Consider the homogeneous system of linear equations formed by the equation yo = 0 and by 
the equations corresponding to the binding inequalities in (llla| - (|llb|) with r = CexpiKp). The 
matrix of this system is formed by the row vectors Uj for j G X U {0} and bj for jej. 

The rank of this matrix is at most p, since y is a non-zero solution to the corresponding 
homogeneous system with p -|-1 variables. Suppose now that the rank of the matrix is less than 
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p. Then, by the Rouche-Capelli theorem (Proposition I3.4L we have a y solution to the system 
such that y and y are linearly independent, and y is an optimal solution to the LP problem (IT^ , 
too. Now, for appropriate real numbers c and c, the vector c • y + c • y is a non-zero optimal 
solution to the LP problem (IT^ , and it satisfies more inequalities in (IllaD - dllbb with equality 
than y does. We can check that the corresponding polynomial Q(:r) := • y^ + c • ym)x"^ 

satisfies all the conditions of the lemma, and it satisfies more inequalities in (I14a|l - (ll4b|l with 
equality than q does, which contradicts our choice of q. Hence the rank of the matrix is precisely 
V- 

Consider again the linear system for the possibly non-zero coefficients arid 

obtained from the equality constraints of the LP problem (IL^ . The matrix of the system is formed 
by the column vectors aj for / € X and bj for j G J", and the right hand side of the system is just 

aj. Since this system is consistent, the Rouche-Capelli theorem implies that the rank of its matrix 
is equal to the rank of its augmented matrix. We already know that the rank of its augmented 
matrix is precisely p, from which the proof of the statement follows. □ 

Lemma 3.11. Let positive integers k and p > 1 be such that Cinip{k,p) > 0. Assume that and 

are the coefficients of an implicit LM method of order p with k steps and with SSP coefficient 
Cimp(k, p), and let q he any non-zero polynomial of degree at most p that satisfies the system (I14all - (ll4d|) 
with r = Cimp{k,p). Let us define the sets 1 ■={]■. 1 < j < k, and q{j) = 0} and LT 0 < 

i < k, and — q'd) Cimp{k, p)q{i) = 0}, and let us introduce the vectors aj := {i'")m=of‘^r 0 < j <k, 

and bj := ^imp{k, p)f^)^^Q for 1 < j < k, and bo = (0, —1,0,.. .,0). Then the following 

statements hold. 

(i) If 5j f Ofor some 1 < j <k, then q{i) = 0. 

(ii) If f>j 7 ^ 0/or some 0 < j <k, then —q'ij) Cimp{k, p)q{i) = 0. 

(Hi) If the vectors ajforj G 1 and bjforj G LT are linearly independent then there exists no other implicit 
LM method of order p with k steps and SSP coefficient Cj„ip{k, p). 

(iv) The polynomial q can he chosen so that the vectors ajfor ; G X and bj for j E LT span a p dimensional 
subspace o/RP+^. 

Proof. The proof is essentially the same as that of Lemma B.lOl □ 

4 Upper bounds on SSP coefficients 

In this section we derive upper boimds on SSP coefficients using Lemma 13.61 and Lemma 13.71 
Proposition 14.11 is just the classical upper bound on the SSP coefficient for explicit LM methods, 
found in IILen89l . 
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Proposition 4.1. Let k, p be positive integers. Then the following inequality holds. 

ifk>2 and p <k 
Cexp{k/ p} ^ 0 if p ^ k 

U ifk = p = l. 

Proof. Let us consider q{x) := x{k — x)P^^, a polynomial of degree p. If ?c > 2 then q satisfies the 
system (I14all - (ll4cll with r = max If ^ = 1 and p > 2 then q satisfies the same system 

with r = 0. Finally, if p = A: = 1 then q satisfies the same system with r = 1. The statement of the 
proposition follows from Lemma 13.61 □ 

Theorem 14.11 can be viewed as a refinement of the upper bound Cimp(A, p) < 2 for all k > 
1 and p > 2 proved in IILen91llHR05l . 

Theorem 4.1. Let k>l and p >2be arbitrary integers. Then the following inequality holds. 

ifk > 2 and p <2k 
if p > 2k 

ifk = l and p = 2. 

Proof. Let us consider q{x) := x^{k — xY~^, a polynomial of degree p. If A > 2 then q satisfies the 
system (I14all - (ll4db with r = max . If A = 1 and p > 3 then q satisfies the same system 

with r = 0. Finally, if A = 1 and p = 2 then q satisfies the same system with r = 2. The statement 
of the theorem follows from Lemma 13.71 □ 


f-imp {k/ p) ^ 


2k—p 

TT 

0 

2 


Proposition 4.2. The following statements hold. 

(i) Suppose that positive integers p and A are such that Cexp(A, p) > 0. Then there exists an explicit LM 
method of order p with A steps, and with C = Cexp(A, p). 

(ii) Suppose that integers p >2 and k > 1 are such that Cimp(A, p) > 0. Then there exists an implicit 
LM method of order p with A steps, and with C — Ciinp(A, p). 

Proof. The proof of statement (HJ is based on Lemma l3Tl and Proposition 14.11 and on the fact that 
if a polynomial q solves the system (IT^ with r = ro > 0, then the polynomial q := q — q{0)/2 
solves the same system with r = ri for some 0 < ro < ri. Similarly, the proof of statement (HU is 
based on Lemma |T2] and Theorem 14.II □ 

Theorem 14.21 presents a relation between the optimal SSP coefficients for explicit and implicit 
LM methods of different orders. This inequality, together with the existence of arbitrary order 
implicit LM methods with C > 0 proved in IISan86L implies the existence of arbitrary order LM 
methods with C > 0. Using this approach, we can show that Cexp(A, p) > 0 for a given order p, 
if the number of steps is at least A > 2P^^ + 1. In contrast. Theorem 15.11 proves the existence of 
explicit LM methods with C > 0 for a much lower number of steps. 
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Theorem 4.2. Let k, p be positive integers. Then the following inequality holds. 

CimpiKlp) <2Ce^p{Kp). (17) 

Proof. Lemma 13.61 guarantees the existence of a non-zero real polynomial q of degree at most p 
that satisfies the system (I14a|l - (ll4cb with r = C^xpiK p). Let us consider the polynomial q := q^. 
Clearly q is a non-zero real polynomial of degree at most 2p. We can check that q satisfies 
conditions (|14a|l . (I14cb and (|14d|l . 

The above q satisfies condition (Il4bb . too, with r = 2Cexp(k/P)- This is because —q'ij) + 
2Cexpik,p)qii) = 2q(i) {—q'ii) + CexpiKp)q{j)) > 0 for all 1 < j < k. Here we have used that q 
satisfies conditions (I14ab and (Il4bb with r = CgxpiKp)- Hence q satisfies the system (I14ab - (ll4db 
with r = 2Cexp(fc, p), and thus Lemma ISTI implies the inequality (fT7)l . □ 


5 Existence of arbitrary order SSP explicit LMMs 


In IILen89[ Theorem 2.3(ii)], it is asserted that there exist explicit contractive linear multistep meth¬ 
ods of arbitrarily high order; the justification cited is IISan86l Theorem 2.3]. The latter Theorem 
does prove existence of arbitrary-order contractive linear multistep methods; however, it uses an 
assumption that if ocj f 0 for some /, then also fj f 0, which cannot hold for explicit methods 
(since necessarily ocq 0 and fo = 0). Hence the methods constructed there are necessarily 
implicit. 

In view of Lemma 13.61 the existence of arbitrary-order SSP explicit LMMs can be shown by 
proving the infeasibility of (|14a|l . (Il4b|l and (I14c|) with r = 0 and a large enough k for all p. 
This can be achieved by applying Markov brothers-type inequalities to the polynomial q. These 
inequalities give boimds on the maximum of the derivatives of a polynomial over an interval in 
terms of the maximum of the polynomial, and they are widely applied in approximation theory. 

Proposition 5.1 (Markov brothers' inequality IIMar90 1 lMar92 II 1. Let P be a real univariate polynomial 
of degree n. Then 


max < 


n^{n^ 


1^) (n^ — 2^) • • • (n^ — — 1)^) 

1-3-S'--(2^-1) 


max \P(x)\ 

-1<X<1 


for all i positive integers. 

Theorem 5.1. Let positive integers k and p be such that 


k> 


/p2(p2_l) 

p + 2 

/ 6 

2 


(18) 


Then Cexp(k, p) > 0, that is, there exists an explicit LM method with k steps, order of accuracy p and 
C>0. 
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Proof. Suppose to the contrary that Cexp{k,p) = 0 for some p and k satisfying dlSl) . Then by 
Lemma |3]6j there exists a non-zero real polynomial q of degree at most p that satisfies the system 
dMaMMcD with r = 0. We choose this q to satisfy properties [T] and |6] of Lemma 13.81 That is, q is 
of degree p, and q{x) > 0 for all x G [0, k]. 

First we define a := maX;|.g[o,A:] and b := max^g^jt] \q"{x)\, and we introduce the polyno¬ 
mial P{x) := q{{x -h 1)|) — By construction, max^g[_;^ ;^] |P(x)| = and maXj.g[_;j |P"(x)| = 

^b, and thus the Markov brothers' inequality for the second derivative (Proposition 15.11 with 
£ = 2) applied to P implies that 


k^ p2(p2-l) 

b ■ — < a - ^ 


(19) 


Now we give an upper estimate for a in terms of b. By the definition of b, and since qf]) > 0 
for all integers l<i <k due to condition (I14b|) . we have q'{x) < b{l — x) for x G [0,1], and 
q'{x) < b ■ max(x — j,j -|-1 — x) for x G [/',/ -|- 1] for all integers i<;<fc-i. From this, we get 
that Jg q'(x) dx < I and dx < | for all for all integers 1 < J < k — 1. 

Since q' is a polynomial of degree p — 1, it follows that 


G Z : 0 < ; < fc 


1, and 


max q'{x) > 0 > 
^eTi+i] J 


< 


I 

12 


that is, q' can take positive values on at most subintervals of [0,k] of length one. Since 
q{0) = 0 and q is non-negative on [0,k], the interval [0,1] is one of these subintervals. 

Now, because of <^(0) = 0, we can estimate a from above by 


a < 



max((j'(x),0) dx < - -|- 


b 

4 


p-2 

b 

P + 2 

2 j 

4 

2 


Combining and (ll^ yields that 


( 20 ) 


k < 


/p2(p2-l) 

p + 2 

/ 6 

2 


which contradicts the assumption of the theorem. Hence CgxpiKp) > 0 for all k and p positive 
integers satisfying dlSll . □ 


6 Asymptotic behavior of the optimal step size coefficient for large 
number of steps 

Since for all positive p, C^xpiKp) and Cimp(^/P) are trivially non-decreasing in k, the following 
definitions are meaningful. 

Cexp(®®/P) ■— lirri Cexp(^/p) f'imp(®®/P) •— Hm Cimp(^/p) 

k^co ^ ^ Ic-S-oo 
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Due to the existence of arbitrary order explicit and implicit LM methods with C > 0 (Theorem 
I5.1D , and due to the upper bounds on SSP coefficients (Proposition 14.11 and Theorem 14.IL 0 < 
Cimp(oO/ p) < 1 for all positive p, and 0 < Ciinp(oo, p) <2 for all p >2. 

Theorem 6.1. Let p be a positive integer. Then the following statements hold. 

(i) If p is odd then there exists a positive integer Kp such that for all integers k> Kp 


Cexp{k/ p) — Cexp(®®/P)- 


( 21 ) 


(ii) If p is even then there exists a positive integer kp such that for all integers k > kp 

^ _ 2 — f'exp(oO/ P ~ 1) ~ Cf.xp{k, p) < p • 

Proof. For positive integers p and n < p, let L{p,n) denote the smallest whole number L such 
that for all k large enough that C^xpiK p) > 0, there exists a polynomial q{x) := c • Ylm=o(^ ~ 
of degree p that satisfies the system (ll4all - (ll4cD with r = Cexp{k, p), and whose roots satisfy 

\{m\l <m< p,Re(A„) G [0, L]}| > n; (23) 

if no such L exist, then we define L{p,n) := oo. Let N{p) denote the largest n < p positive integer 
such that L{n,p) < oo. Because of properties [THl] and |4] of Lemma |3]8j N(p) is odd for all p > 1. 
We break the proof into a number of steps. 

Step 1. Let integers k >2 and p > 1 be such that Cgxp{k, p) > 0, and let cj he a polynomial 
that satisfies the system (I14all - (ll4cll with r = Cexp{k, p). We can check that q{x) := {k — x) ■ q{x), 
a polynomial of degree p + 1, satisfies the same system with r = max (O, ro — ^), and thus 
Lemma 13.61 implies that CgxpiKp + 1) < max (O,Cexp(fc, p) — i^)- Hence Cexp(k, p) is strictly 
decreasing in p until it reaches 0. 

Step 2. Suppose that N(p) = p for some positive integer p. We choose the positive integer Kp 
so that CexpiKp,p) > 0 and Kp > L{p,p) + p/Cexp{Kp,p). Such a Kp exists, because Cexp{k,p) is 
non-decreasing in k. Let qhe a polynomial of degree p that satisfies the system (IMaMMcIl with 
f = Cexp{Kp, p) and k = Kp, and whose roots satisfy (1^ with L = L{p, p) and n = p. We can 
check that q satisfies the inequality (I14all for all j > Kp. In view of statement (|ii]l of Lemma 13.31 
and because of 


P 1 

m=l J 


P I 


< p 


i-L{p,p) 


^ f'exp(fxp/ p)/ 


q satisfies the inequality (Il4b|l with r = Cexp{Kp, p) for all j > Kp. 

Therefore q satisfies the system (I14all - (ll4cll with r = Cexp{Rp,p) for all k > Kp, and thus 
Lemma 1331 implies that Cexp(k, p) < Cexp{Kp, p) for all k > Kp. Since Cexp{k, p) is non-decreasing 


19 



























in k, this means that C^xpiKp) = Cexp{Kp,p) for all k > Kp. Hence for all positive p, N{p) = p 
implies that condition (|2T1) holds for an appropriate choice of Kp. 

Step 3. Suppose that N(p) = p holds for some positive integer p, and suppose that q is a 
polynomial of degree p that satisfies the system (ll4aD - (ll4cll with r = Cexp{^,p) for all positive 
k, and whose roots satisfy (|2^ with L = L{p,p) and n = p. Such a q exists in view of Step 
2. Let us define p := p/Cexp(oo,p), and let us introduce q{x) := (x^ — 2 {^p + L(p,p)) • x + 
(4p + L{p, p))^ +4p^) • q{x), a polynomial of degree p + 2. 

We can check that q satisfies the inequality (I14all for all positive j, and that q satisfies condition 
(ll4cD . Using the inequality !/(/ — A) + !/(/' — A) < 1/ Im(A), we can show, in a similar way as 
in Step 2, that q satisfies the inequality (I14bl) with r = |Cexp(oo, p) for all j > L{p, p) + 2p. Using 
that 1/(/ — A) + 1 /(j — A) is decreasing in j in the range 0 < j < Re(A) — Im(A), we can show 
that q satisfies the inequality (|14b|l with r = Cexp(oO/P) — 2(L(p, p) +4p)/ {(L{p,p) +4p)^ +4p^) 
for all 1 < i < L{p,p) + 2p. Therefore q satisfies the system (ll4all - (ll4cD with the maximum of the 
above expressions for r for all positive k, and thus Lemma 13.61 implies that 


Cexp(oo, p + 2) < max 


r^exp (oo,p), Cexp (oo,p) 


2(L(p,p)+4p) \ 

(L(p,p) +4p)2 +4p2y 


Hence for all positive p, N(p) = p implies that Cexp(oo, p) > Cexp('^/ P + 2). 

Step 4. Suppose that N(p) = N(po) = po for some p > po positive integers, and let e > 0 be 
arbitrary In view of Step 2, there exists a positive integer Kpg such that CgxpiK po) = Cexp(oO/ Po) 
for all k > Kpg. Let us define L := max(L(p, po), Kpg) and fC := L + (p — po)/£. 

Since N(p) = po, we can choose k > K so that there exists a q{x) := c ■ ~ 

polynomial that satisfies the system (I14all - (ll4cll with r = Cexp(k, p), and whose roots satisfy 
the conditions Re(Ai),...,Re(Ap(,) G [0, L] and Re(App+i),... ,Re(Ap) > K. Let us introduce 
q := c ■ (—l)P^Po . “ ^m), a polynomial of degree po- We can check that q satisfies the 

conditions of statement (E]) of Lemma 13.31 for all positive j, and that q satisfies condition (I14cll . 
Moreover, q satisfies the conditions of statement (El) of Lemma 1331 with r = Cexp(k/ p) + £ for all 
1 < j < L, since 


Pi Po 1 Pi Po -I 1 Po 1 

m=0/ m=of m^po+lJ m=0 J ^ ^ m=01 

for all 1 < / < L. Now Lemma 13.31 implies that q satisfies the system (ll4aD - (ll4cll with r = 
Cexp(k, p) T £. From this. Lemma [3^ implies that Cexp(k, p) T £ ^ Cexp (P^Po) = Cexp {Kpg,p0) = 
Cexp(oO/Po)- Since e > 0 is arbitrary and Cexp(oO/P) is non-increasing in p, this means that 
Cexp{oo,p) = Cexp(oo,po). Hence N(p) = N(po) = po implies that Cexp(oo,p) = Cexp(oo,po)- 

Furthermore if p = po -|- 1 then p is even, and thus property [5] of Lemma guarantees that 
Ap = k. Then a similar reasoning yields that N(p) = N(p — 1) = p — 1 implies Cexp (oo, p — 1) — 
CexpiK p) < -^ for all k> L. 

The proof of statement (E)) of the Theorem. Clearly, N(l) = 1 holds, due to condition dlTcIl . 
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Assume now that N{2j + 1) =2] + 1 holds for all integers 0 < / < m — 1. Since 1 < N{p) < p 
and N{p) is odd, we have N {2m + 1) = N(2/o + 1) for some integer 0 < jo < m. Suppose that 
jo < m. Then Step 4 implies that Cexp(2?n + 1) = Cexp(2/o + 1)/ which contradicts the result of 
Step 3. Therefore N{2m + 1) = 2m + 1. 

By induction we conclude that N{p) = p for all odd positive integers p. Finally, Step 2 implies 
statement (0 of the Theorem. 

The proof of statement © of the Theorem. Let p be an arbitrary even positive integer. Then 
N{p) = po for some odd positive integer 1 < po < p, since 1 < N{p) < p and N{p) is odd. We 
already know that N{po) = po, thus Step 4 implies that Cexp(oo, p) = Cexp(oO/Po)- Suppose that 
Po 7^ P “ 1- Then Step 3 implies that Cexp(oO/ p) = Cexp(oO/ Po) > Cexp(oO/ P — 1)/ which contradicts 
the fact that Cexp(oO/ p) is non-increasing in p. Therefore po = p — 1, and thus Step 4 implies that 
Cexp (oo/ P “ 1) “ Cexp{k, p) < for a fixed kp and for all k > kp. 

For k large enough, the inequality Cexp (oo/ P “ 1) “ Cexp(^/ p) follows from Step 1 and 


the fact that Cexp(k/ P — 1) is non-decreasing in k. □ 

Theorem 6.2. Let p >2be an arbitrary integer. Then the following statements hold. 

(i) If p is even then there exists a positive integer Kp such that for all integers k> Kp 

Cimp{k,p) = Cimp P)- (24) 

(ii) If p is odd then there exists a positive integer kp such that for all integers k > kp 

jfZTl — P ~ 1) ~ Cjmpik, p) < ■ (25) 

Proof. The proof is essentially the same as that of Theorem 16.II □ 


The proofs of the following two theorems make use of the assumption that step sizes are 
fixed. 

Theorem 6.3. Let p be an odd positive number, and let Kp be the smallest positive integer such that 
Cexp{Kp,p) = Cexp(oO/P)- Then for all k > Kp there is a unique optimal explicit LM method with k 
steps, order of consistency p and SSP coefficient Cexp(oo, p), and for different values of k, the non-zero 
coefficients of the unique optimal methods coincide. Moreover, there is a finite sequence of positive integers 

jifji, ■ ■ swc/z that 1 < /I < ;i -|- 1 < ;2 < ;2 -l- 1 < ... < jvzl < izzl + 1' coefficients 

2 2 2 
of the optimal method satisfy 5j = 0 for all 1 < j < k, and that if f>j 0 for some 1 < j < k then 

i £ {1/ ji ,/i + 1/ • • • / jEM'lEf: + !}■ 

Proof. Suppose that p and Kp satisfy the conditions of the theorem. Theorem 16.II guarantees the 
existence of such a Kp for all odd positive integers p. Let us set Kp := 2Kp. Then Cgxp{Kp,p) = 
Cexp(oO/ p) > 0, and in view of Lemma |3A] and property [T] of Lemma lT8l there exists a q polyno¬ 
mial of degree p that satisfies the system (ll4a|) - (ll4cD with r = Cgy^pjoo, p) and k = Kp. 


21 









First we show that q satisfies the inequality ([Ma) for all 1 < ;■ < Kp with strict inequality. 
Suppose to the contrary that q{jo) = 0 for some positive integer 1 < jo < Kp. If jo > Kp then q 
satisfies the system (I14all - (ll4cll with r = Cexp(oo, p) = Cexp{Kp, p) and k = Kp, but does not satisfy 
property |2] of Lemma 13.81 which is a contradiction. Similarly, if jo < Kp then the polynomial 
q{x) := q{x — jo) satisfies the system (I14all - (ll4cll with r = Cexp(oo,p) = Cexp{f^p,p) and k = Kp, 
but fails to satisfy property |2] of Lemma 13.81 which is again a contradiction. Hence none of the 
inequalities in (I14all hold with equality. 

Since —q' + Cexp(oo, p)q is a polynomial of degree p, it satisfies at most p inequalities in (|14b|l 
with equality. Because of statement dr^ of Lemma 13.101 and due to q(i) >0 for all 1 < ; < 
Kp, we can choose q so that the number of binding inequalities in (Il4b|l is precisely p. Then 
—q'{0) + Cexp(oO/p)(?(0) = —q'{0) < 0 because of properties |4| and |6] of Lemma 13.81 Let J' and 
bj for 1 < j <Kp be as in Lemma 13.101 We can see that (I14bl) holds for all 1 < / < Kp if 

and only if J = + 1,.. ■ fjEzl'jpzl +1} a sequence of integers ji,j 2 , ■ ■ ■ ,jp-i satisfying 

2 2 _ “ 2 “ 

1 < /i < /i + 1 < ... < jp -1 < jp -1 +1 < Kp. Now property |2] of Lemma |3^ implies that q{j) > 0 

for all j > Kp. We can see that —q'{j) + Lexp(oo, p)^(/) > 0 for ] > Kp. This is because the set 
J' contains all the p roots of —q' + Cexp(oO/ p)q- It follows that q satisfies the system (I14all - (ll4cll 
with r = Cexp(oO/ p) for all positive k. 

Suppose that 5j and fij for 1 < j < k are the coefficients of an explicit LM method with 
k > Kp steps, order of consistency p, and SSP coefficient Cexp(oO/P)- Proposition 14.21 guarantees 
the existence of such a method for all k > Kp. Statement (0 of Lemma 13.101 implies that Sj = 0 
for all 1 < / < fc, and statement dull of Lemma 13.101 implies that if fij = 0 for some f < j < k, 
then j G J. Since bj for j ^ J are linearly independent (by statement (|i0 of Lemma 13.101 and 
our choice of q), statement dHIll of Lemma 13.101 implies that there is a unique optimal method 
with SSP coefficient Cexp(oO/ p) for all k > Kp. The possibly non-zero coefficients of this unique 
optimal method, namely jij for j G J7, satisfy the same linear system for all k > Kp. Thus for 
different values of k, the non-zero coefficients of the unique optimal methods coincide. □ 

Theorem 6.4. Let p be an even positive number, and let Kp be the smallest positive integer such that 
Cimp(f^p/P) = Cimp{oo,p)- Then for all k > Kp there is a unique optimal implicit LM method with k 
steps, order of consistency p and SSP coefficient Ciinp(oO/P), and for different values of k, the non-zero 
coefficients of the unique optimal methods coincide. Moreover, there is a finite sequence of positive integers 

jirji, ■ ■ ■ ,jvzir such that l<ji<ji + l<j 2 <j 2 + 'f<---< jp -2 < jpzl + 1/ that the coefficients 
2 “ 2 “ 2 
of the optimal method satisfy Sj = 0 for all 1 < i < k, and that if Bj 0 for some 0 < i < k then 

j e{0,l,ji,jr + l,...,j^,j^+l}. 

2 2 

Proof. The proof is essentially the same as for Theorem l6.3l □ 
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